RAL-TR-97-039 



What Have We Learned About and From Fl{x,Q^) at HERA? 



Robert S. Thorne 

Rutherford Appleton Laboratory, 
Chilton, Didcot, Oxon., 0X11 OQX, U.K. 



Abstract 

Recently the HI conaboration has pubhshed a "determination" of the structure function 
Fl{x, Q'^) at low X. I address the question of how reliable this determination really is. I argue that 
it is in fact a consistency check of a given theoretical approach rather than a real determination 
of the value of Fl{x, Q^), but potentially a very useful one. I compare the consistency of different 
approaches, and indeed find that a LO-in-Os calculation of structure functions is completely ruled 
out. I also find that the "determined" values of Fi^{x, Q^) are surprisingly stable under changes in 
theoretical approach but, when working consistently within a well-defined theoretical framework, 
the values of Fl{x, Q^) implied are somewhat lower than previously quoted. 
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1. Introduction 

The recent measurements of F2{x,Q'^) at HERA [1][2] have provided data on a structure 
function at far lower values of x than any previous experiments. However, these measurements are 
not of the structure function directly. Rather they are of the differential cross-section 

{2{1 -y) + y')F2{x,Q^)-y^FL{x,Q^) . (1.1) 

is the squared four-momentum transferred in the lepton-proton scattering, x is the scaling 
variable and y = Q'^/sx where s is the centre of mass energy squared. At HERA s = 90200GeV^. 
However, Fl{x, Q^) is expected to be much smaller than F2{x, Q^) for most of the parameter space 
in which measurements take place (it must be smaller than F2{x, Q^) at all x and Q^), and within 
this parameter space y is nearly always small (0.25 or less). Hence, the measurement is usually 
effectively of F2{x,Q'^). 

However, at the fringe of parameter space, i.e. the highest values for the lowest x mea- 
surements, the values of y can be somewhat larger, reaching pa 0.7 at their maximum. For values 
of y larger than about 0.25 the fact that Fl{x,Q'^) is exprected to have a non-zero value starts 
affecting the value of F2{x, Q^) extracted from the measurement of the differential cross section by 
about 1%. This increases to about 10 — 20% at the highest y values of 0.7 (clearly depending on the 
value of Fl{x, Q^) used). Hence in this region of parameter space, within the typical quoted errors, 
the HERA collaborations cannot be claiming to measure F2{x,Q'^) directly. It would perhaps be 
simplest to release data in this region of parameter space in the form of the differential cross-section 
rather than in terms of structure functions. However, both HI and ZEUS choose to take the values 
of Fl{x,Q^) from some theoretical prediction, and extract the consequent "measured" values of 
F2{x,Q'^). Both collaborations obtain their values of Fl{x,Q^) from fits to the data using the 
usual approach of solving renormalization group equations for parton distributions and combining 
with coefficient functions for parton scattering with the calculations done to NLO-in-as. In [1] 
Fl{x,Q'^) is calculated using the GRV parameterization of partons [3], while in both of [2] it is 
calculated using the parton distributions obtained from their own iterative fit to F2{x,Q'^). 

Thus, both collaborations produce values of F2{x,Q'^) which take into account the values of 
Fl{x,Q^) in an approximate manner. However, it is very important to remember that the values 
of F2{x,Q'^) quoted in [1] and [2] are not really measurements of F2{x,Q^) in the high y region, 
and any attempt to fit the HERA data using some theoretical model should take this into account. 
The values of F2{x, Q^) used in a fit should really be those obtained from the measured differential 
cross-section and the prediction for Fl{x, Q^) consistent with the particular theoretical approach (it 
is difficult to imagine a genuine theoretical approach which does not calculate Fl{x, Q^) alongside 
F2(x,Q^)). Alternatively, the fit could be directly to the cross-section using the calculated values 
of both F2{x, Q^) and Fl{x, Q^). (The two approaches are identical if the error on F2{x^ Q"^) in the 
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former does not change with the predicted Fl{x, Q^), but is equal to its value when Fl{x, Q"^) = 0.) 
For any "conventional" NLO-in-ag approach the differences will not be all that large, but when 
using different approaches, such as those involving summation of leading logs in (l/x), or even a 
LO-in-CKs calculation, the differences in the predictions for Fl{x, Q^) can be large enough that the 
extracted values of F2{x, Q^) move by amounts greater than their errors on the fringes of parameter 
space. Moreover, the direction of movement tends to be the same for every point, so ignoring this 
effect can have a sizeable effect on the fit. Unfortunately the consistent manner of comparing theory 
to F2{x,Q^) values described above seems not to be done in practically all fits to data. 

Recently the HI collaboration took a rather different approach to the treatment of their high 
y data [4]. Rather than adopt the procedure outlined above, they fitted the data for F2{x,Q'^) 
in the region of low y using NLO-in-CKg QCD, used this fit to extrapolate F2{x,Q'^) into the 
high y region, and from the difference between this extrapolation and the measured cross-section 
obtained a "measurement" of Fl{x, Q^). This determination of the value of Fl{x, Q"^) relies only on 
assuming the correctness of the NLO-in-ag fit for F2{x,Q'^). But if one is to assume this correctness 
then it seems perverse indeed not to assume the correctness of the NLO-in-Og Fl{x,Q'^) that is 
predicted as a result. Hence, if the theoretical approach is correct, Fl{x,Q^) is already determined. 
The difference between the measured cross-section and the extrapolated F2{x,Q'^) then provides 
nothing more than a consistency check: if it disagrees with the predicted Fl{x,Q'^) it suggests 
that the theory is wrong, if it agrees then the theory is not necessarily wrong. This is all that one 
can conclude. There may well be other theoretical approaches which fit the low y data equally 
well (or better) but which have different extrapolations and/or predictions for Fl{x,Q^). Indeed, 
in the small x region this is likely since higher orders in are accompanied by higher order 
terms in ln(l/x), and these leading ln(l/x) terms may well introduce important corrections to the 
standard approach. If the extrapolations and predictions match for these other approaches, then 
the predicted Fl{x, Q"^) is just as likely to be correct, but is no more a real measurement than any 
other matching case. 

Hence, the procedure adopted by the HI collaboration is not in itself a real measurement 
of Fl(x,Q^). However, it is potentially a very useful way of discriminating between different 
theoretical approaches. The consistency of the predicted Fl{x,Q^) and the ""measured" Fl{x,Q^) 
for a given approach is a very non-trivial check on the theory since, if the free parameters in 
the theory are already tied down by a fit for high y data, then it is a parameter free check on a 
particular relationship between F2{x,Q'^) and Fl(x,Q'^). Hence, in this paper I will present the 
approach used by HI in the manner of a consistency check. I will do this for four different types 
of approach: a fit analogous to that of HI in [4], which I will demonstrate is not actually a useful 
NLO-in as fit: a more correct NLO-in-ag fit: a LO-in-a^ fit: and a fit using the LORSC approach 
which includes leading \n{l/x) terms[5]. I will also include more data at high y than that in [4]. I 
will find that the "measured" Fl{x,Q'^) is very similar in each of the three latter approaches, but 



2 



lower than in the first one. However, this similarity would start to disappear at higher y. I also 
find that find that the data with its current errors is not yet sufficient to rule out any approach 
other than the LO-in-CKg fit, which it does rule out very clearly. One can, if one wants, interpret 
the "measurements" of Q"^) within the two approaches satisfying consistency to be the likely 

values of FLixjQ"^), but this is not really fundamentally different from quoting instead the values 
of Fl{x,Q^) predicted within either of these approaches. 



2. Comparison of Different Theoretical Approaches. 

I compare the four different approaches outlined above. First I discuss the data to be used 
and the accuracy of Fl{x,Q^) achieved. I define a rescaled differential cross-section 



^ ^ " 27ra2 (2(1 - y) + y^) dxdQ^ ■ ^ ^ 



Hence I may write 



a{x,Q') = F,{x,Q')-Y-'.FUx,Q'), y=ML^±^. (2.2) 

This is clearly a useful definition since in the limit Fl{x,Q^) or y ^ 0, 5- = F2{x,Q^). 
Hence it allows a simple parameterization of the effect of non-zero Fl{x,Q'^) on the extracted 
F2{x,Q^). From (2.2) it is clear that denoting F2{x,Q'^) obtained from the the fit by i^l* (x,(5^), 
the "measured" value of Fl{x,Q^), denoted by F]^{x,Q'^), is given by 

F^ix, Q2) = y(F|(x, Q2) _ aix, Q^)), (2.3) 

and that the error in Fj^{x, Q^) is 

1 

AFr(x,Q2) = y- [(Aa(x,Q2))2 + (AF|(x,Q2))2]2^ (2.4) 

where Aa{x,Q'^) >> AF2{x,Q'^) in practice. Hence, the error on F^{x,Q^) depends linearly on 
the error on the measurement of the cross-section and roughly quadratically on y~^. Although 
this means that, if one wishes to make a sensible consistency check for F^{x,Q^) (by comparing 
to the predicted F^{x,Q^)), it is necessary to have high y, a relatively small decrease in y can 
be countered by a large decrease in the error on the measurement of a{x,Q'^). For example, an 
error of 5% on the measurement of a{x, Q"^) at y = 0.6 leads to a more accurate determination of 
F]^{x,Q'^) than an error of 10% on the measurement of a{x,Q'^) at y = 0.7. Examining the full 
range of data in [1] and [2] there are a number of points which give an accuracy of measurement 
of F]i^{x,Q'^) which is comparable to the 6 points in [4]. Talking a cut on the error produced for 
Fj^{x,Q'^) to be 0.3 and a lower cut on of 5GeV^ (in order not to be too close to the charm 
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threshold), I find that there is one point in [1], and 8 points in the latter of [2]. Of these additional 
points, 5 of the ZEUS points have y comparable to the 0.7 in [4], while the other 3 ZEUS points 
and the HI point have somewhat smaller y but better accuracy in measurement. The full set of 
data points to be used is shown in table 1. 

Now that the relevant data are defined, let me consider the fits. First I note that the heavy 
quark thresholds are treated rather differently than in [4] . Rather than using the NLO-in-a^ fixed 
flavour scheme, where the charm contribution to the structure function is entirely generated by 
coefficient functions convoluted with light parton distributions, I simply change the number of 
active quark flavours discontinuously at m\ = Q^, and treat the heavy quarks as massless above 
this. I take = 2.75GeV^ in order to give a good description of the charm data[6][7], and choose 
m| = 20GeV^. Neither of the above approaches for the treatment of charm is entirely satisfactory 
since the former does not sum leading logarithms in Q^/m^ while the latter does not treat the 
threshold in a correct manner. In practice the former is more accurate for < lOGeV^ and the 
latter more accurate for Q > 20GeV . For lOGeV < Q < 20GeV the approaches are of roughly 
equal validity. 

In comparing the predicted and measured Fl(x,Q^) I adopt a different procedure from HI 
and assume that the variation in (xjQ^) is entirely due to variations within a given theoretical 
approach fitting a particular set of data. The difference between the various theoretical approaches, 
and/or from using different sets of data will be seen in the four different sets of results. Within a 
given approach the variation in -^^(a;, Q'^) is due only to letting the ^or the fit vary to the extent 
that the quality has a confidence level only a little lower than the absolute best fit. I will discuss 
the results of this later, but mention here that the errors arc very small. Also, since I have as much 
faith in the NLO-in-a^ calculation of Ff^{x, Q^) as of (x, Q^), and hence am only performing a 
consistency check, I include all data other than those points in table 1 in the fit, and let the value of 
(5^) used in the extraction of F^{x, Q^) for these points be that predicted by the parameters 
determined in the fit (i.e. the fit is iterative). Thus, the procedure differs from [4] in including the 
small number of points with y > 0.35 not in table 1, which has an extremely small effect on the fit, 
and also by using a particular Q^) for the points in the fit rather than letting it vary over the 

rather extreme range of — > F2{x,Q'^). Avoiding this second variation means that I do not have 
the large systematic error in the extrapolated (x,(5^) which is seen in [4]. Overall the method 
of defining the uncertainties due to the fitting assumptions reduces them immensely compared to 
those in [4]. 

1. First I consider the analogous fit to that performed by HI. In this the fit is performed to BCDMS 
[8] and HI data only. This leads to a problem often encountered (or rather ignored) by those 
attempting to fit small x structure function data: the parameters defining the parton distributions 
which are determined by the best fit turn out to be completely incompatible with some perfectly 
respectable data which have not been included in the fit. The parameters determined by the the 
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fit are not included in [4]. However, performing a fit to precisely the same data I obtained a gluon 
distribution which is far too small at a; > 0.1 to be consistent with the WA70 prompt photon data 
[9], and is even much smaller than the gluon produced by the HI collaboration in the fit in [1]. 
Although there is a certain amount of uncertainty about the true accuracy of the prompt photon 
data, the gluon obtained is very much smaller than any possible lower limit on the data. It is also 
far too small to be consistent with the EMC charm data [6] using any sensible treatment of charm, 
and I note that the gluon required by this charm data is actually similar to that required by the 
WA70 prompt photon data, i.e. ^ 2.5(1 - xf at = 5GeV^ Since there is so little gluon at 
high X in my fit, the momentum sum rule allows much more at small x. This enables the HI data 
to be fit well even using a low value of Aqcd which is preferred by the BCDMS data. Indeed, my 
fit chooses — 186MeV which is very similar to the HI Ag-^ = 210MeV. This should be 

compared with those analyses of only HI data which produce ^qcd ~ 350MeV when the high x 
parton distributions are constrained differently [10]. 

Ignoring this very serious problem, the fit is performed as in [4], and considers only the same 
points as [4] for measuring Fl{x,Q'^). As already mentioned, the only uncertainty in ^2(0;, Q^) is 
that obtained in allowing the quality of the fit to vary. For this fit to a relatively small number of 
data points the variation in while maintaining a "good" quality fit is not insignificant. 

However, since I am not really taking this fit seriously anyway, but are only using it for illustrative 
purposes, I quote the error as being due to the cross-section measurement only. As one would 
expect, the results of this procedure are very similar to those in [4]. The fit to the small x HI data 
is very similar, as is the extrapolation. The values of FJ^{x, Q^) for the 6 data points for a used in 
[4] are shown in table 2 along with the values "measured" in [4](denoted by F^{x, Q'^)[4\) and the 
values of of F^{x,Q^). The results are also displayed in fig. 1, where they are compared directly 
with those in [4]. One can see that the results for FJ^{x,Q^) are indeed very similar to those in 
[4]. The predicted Fl{x,Q'^) is very slightly higher than in [4] at lowish but the two converge 
at higher Q^, which can be interpreted as the effects of the different treatment of charm in the two 
approaches not compensating exactly at lowish Q"^^ but the difference disappearing at higher Q^, 
as expected. Hence, I have very good compatibility with [4] so far. 

I am now able to explore the implications of this type of fit quantitatively. I do indeed have the 
previously mentioned dramatic inconsistency with the prompt photon data and the EMC charm 
data. Also, if I simply put the parameters obtained by the previous fit into a global fit to structure 
function data, i.e include NMC [11], CCFR [12] and ZEUS data, then I find that the fit to ZEUS 
data is a little worse than that obtained form a global NLO-in-a^ fit, and the fit to NMC and 
CCFR data is very much worse. Overall, the global fit is completely uncompetitive with the 
fits produced by e.g. MRS or CTEQ even ignoring the particular problem of the high x gluon. 
Hence, the type of fit performed above certainly does not lead to a correct parameterization of the 
parton distributions and hence extrapolation of the structure function. Thus, even assuming the 
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correctness of the NLO-in-ag calculation of structure functions one does not really learn anything 
concrete about Fl{x,Q'^) from the above fit, or from [4]. 

2. In order to make a more reliable consistency check for the high y data I now consider the 
implications of using precisely the same theoretical framework but including more data in the 
fit. In order to see the correct results of using the NLO-in-ag calculation I must determine the 
free parameters in the calculation by fitting to all available structure function data, with some 
appropriate cuts, and also imposing the constraint that the charm structure function data from 
EMC [6] and HI [7] is also described well. In practice I perform the fit to data in precisely the 
same manner as in [5] except for one minor point. The ZEUS shifted vertex data in the latter of [2] 
seems to be systematically larger than the data in the former. Hence I allow their normalization 
to be 0.98 that of the nominal vertex data, i.e. the amount allowed by the uncorrelated errors on 
normalization (the fit would actually prefer 0.97). Hence, the fit is almost identical to the NLOl 
fit in [5], but this shift in normalization allows the to be a few points better. It is qualitatively 
a great deal different from the fit discussed above. The constraint on the gluon at high x results in 
there being much less gluon at small x and the value of Ag^T^ rises to 315MeV in order to produce 
a satisfactory fit to small x data. This rise in Ag^^ results in a worse fit to BCDMS data, but the 
total fit quality is a great deal better than that for the type of parameterization above. In the fit 
the normalization required by the HI data is 0.985 and for the ZEUS data it is 0.99. The quality 
of the absolute best fit has a rather low confidence level (assuming Gaussian errors) of 6%. Letting 
this reduce to 1% results in a variation of (x, Q^) which when inserted in (2.4) leads to negligible 
effect and I take the error to be that due to Aa{x, Q^) alone. 

Using this fit the form of -F|'(x, Q'^) is rather different at small x from that obtained using the 
first procedure. The growth with is slower than above, even taking into account the different 
normalization of the data, and hence in a given x bin the value of -^^(x, Q^) becomes progressively 
smaller as y increases. From (2.3) it is obvious that this will result in a lower value of F]^{x,Q'^). 
The results for F^{x,Q'^), along with F]^(x, Q^) 

are shown in table 3 and displayed in fig. 2 where 
now I use all 15 data points. (I display the extracted points with a negative value of F^{x,Q'^) 
if this is what is obtained from the use of (2.2), even though this is forbidden in practice.) In the 
top figure all data with y ~ 0.7 is shown, while the lower figure contains the points with rather 
lower y. The measured values are indeed smaller than previously. Compared with the values in [4] 
F^{x,Q'^) is about 0.1 smaller for the lowest x points increasing to more than 0.15 lower for the 
highest X and points. Hence, for some of these points the value of F][^{x, Q'^) reduces to less than 
60% of its value in [4]. By examination of fig. 2 one sees that, as well as the values of F£^{x, Q^) 
for the 6 points corresponding to those in [4] having reduced significantly, the inclusion of the 
other data points tends also to imply a somewhat lower value of F^{x, Q"^). A simple unweighted 
average of the 6 points in [4] gives F^{x,Q^) = 0.48. The same 6 points in this extraction yield 
F^{x,Q'^) = 0.34, while for all 11 points with y f« 0.7 I get F^{x,Q'^) = 0.27. I claim that if one 
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believes the validity of the NLO-in-ag calculation then the values in table 3 are the most reliable 
extractions of F|"(x,Q^) since they result from a much more realistic fit. However, taking these 
values results in a rather different quantitative conclusion on the value of Fj^{x,Q'^) for y 0.7 
for 0.001 > X > 0.0001 from that reached in [4], i.e. it is «i 60% of that value. 

Also shown in table 3 and fig. 2 is the theoretical prediction Fl{x,Q^) for this fit. As in 
the previous fit, and as in [4], it is calculated using the NLO-in-ag parton distributions obtained 
from the fit to F2^{x,Q'^) and the order and order longitudinal coefficient functions[13]. 
Hence, I use coefficient functions at one higher order for the longitudinal coefficient function than 
for F2{x,Q^) because they begin at one higher order. This is the correct procedure to obtain an 
expression consistent with the use of the two-loop coupling constant in a given renormalization 
scheme. Ff^(x,Q'^) is a little smaller than that in the first approach in the region shown (it is 
actually larger for x > 0.001), with the difference becoming larger at smaller x. This qualitative 
result is due to the gluon being rather flatter at small x than the gluon in the former approach, 
which grows quickly as x decreases. As one can see from fig. 2 there is reasonable agreement 
between the prediction and the measurement, even though there seems to a systematic trend for 
the values of F]^{x, Q'^) to fall as x increases while over most of the range the values of Ff^(x, Q'^) 
increase slowly with increasing x. The comparison of FJ^{x, Q'^) to F^{x, Q'^) yields a of 13.9 for 
15 points. Thus, this consistency cheek gives no real reason to distrust the NLO-in-Os approach. 
However, it may be argued that both the quality of the global fit and the comparison of F^(x, Q^) 
with F^{x, Q^) in the high y region do not give overwhelming support for this approach. 

3. I now consider a fit to precisely the same data but using a different theoretical framework. I use 
the LORSC calculation of structure functions discussed in [5], which includes all term which are of 
lowest order in ag for all different types of term in x, and does so in terms of physical quantities 
rather than unphysieal, definition-dependent parton distributions and coefficient functions (see 
[14] for a definition of physical anomalous dimensions). Hence it includes a summation of the 
so-called leading ln(l/a:) terms for physical quantities. It is only done to leading order because 
the full set of terms required for the NLO calculation is not yet known. The fit uses leading order 
Aq-^^ = lOOMeV, which corresponds to as(M|) = 0.115. In this fit the normalization required 
by the HI data is 1.00 and for the ZEUS data it is 1.015. The quality of the absolute best fit has 
a higher confidence level than that of the NLO-in-ctg fit, i.e. 34%. Letting this reduce to 10% 
results in a variation of -Ff (x, Q^) which again when inserted in (2.4) leads to negligible effect, and 
I take the error to be that due to Aa{x,Q'^) alone. If I were to let the confidence level reduce to 
1%, as in the NLO-in-ttg fit then the uncertainty in F^x, Q^) would start to increase the error on 
F£^{x,Q'^) significantly. 

Using this approach the values of F^{x,Q'^) in the region of x under consideration are about 
1.015 times those in the NLO-in-ag approach at = 6GeV^. Taking into account the difference in 
normalization of the data in the two approaches this means that when using HI data the comparison 
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to data is very similar for the two approaches at this Q^, while for the ZEUS data the LORSC 
F2{x, Q^) is slightly smaller relative to data than the NLO-in-ag prediction (due to the 2.5% shift 
in normalization of the data). The growth of -FfC^j Q^) with is greater in the LORSC approach, 
even taking into account the normalization difference, and this leads to the better global fit, as 
discussed in [5]. This greater growth means that the values of Fj^{x,Q^) extracted using this 
approach start out very similar to the NLO-in-ag approach for HI data and not much above 
6GeV^, but steadily increase with respect to them as increases. For the ZEUS data the larger 
normalization shift leads to the LORSC values of Fl{x,Q'^) being slightly smaller than those for 
the NLO-in-CKg fit for lowish Q^, but they increase steadily above those for the NLO-in-ag fit at 
higher Q'^. The values of F^{x,Q'^) are shown in table 4 and fig. 3 and by comparing to table 
3 and fig. 2 the trend described above can be easily observed. It results in the average value of 
-F™(a;, Q^) being slightly higher in this approach than the previous one, and in the tendency of 
the values of FJ^{x,Q'^) to fall with increasing x or to reduce. The discrepancy between the 
extracted values for the two theoretical approaches, though not large for the points considered in 
this paper, would continue to increase for increasing Q^. 

Also shown in table 4 and fig. 3 is the theoretical prediction F^(x, Q^) for this fit. As discussed 
in [5] it is smaller than that in the NLO in approach in the region shown but the difference 
becomes less at smaller x, where the NLQ- in-ag curve at constant y = 0.7 falls as x goes to 0.0001 
and below, while in the LORSC approach it stays more constant. As one can sec from fig. 3 there 
is again reasonable agreement between the prediction and the measurement, though some of the 
values of F^(x,Q^) at smaller x lie quite a long way above the F^{x,Q^) curve. The comparison 
of Q^) to F^{x, Q^) yields a of 12.4 for 15 points. This is very slightly better than in the 

NLO-in-Os approach. Hence, this consistency check is certainly satisfactory, giving no evidence for 
the failure of the LORSC approach. Indeed, it is even slightly better than the NLO-in-ag approach 
which has a rather less successful global fit to data. 

4. Finally I consider the consequences of using a simple LO-in-Ocj calculation of the structure 
functions. Methods of fitting the small x data based simply on the LO-in-Og calculations have 
not been uncommon in recent years, the "double asymptotic scaling" formula comes from this 
calculation [15], and have been advertised as being relatively successful, especially if is not too 
low. I find that the standard approach is rather too simplistic, and in fact the LO-in-Og approach 
is completely ruled out. If one performs a global fit to structure function data other than that used 
in the extraction of -F™(x, Q^) and using simply those values of F2{x, Q"^) advertised in [1] and [2], 
then the result, requiring the same normalizations as the NLO-in-Os approach, is very nearly as 
good as for the NLO-in-ajg fit over the whole x range. It yields a of ~ 20 more for ~ 1100 
points and a confidence level of 2%, compared to the earlier 6%. However, the quality of the fit 
to the small x data is achieved at the expense of having a fairly large coupling constant, one-loop 
as{M^) = 0.124 and a much larger small x gluon than the NLO-in-a^. This, coupled with the 
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fact the the NLO corrections to Fl{x, Q^) are negative at small x, leads to Ff^{x, Q^) being a great 
deal larger at small x than in the NLO-in-ag approach. Moreover, as in the NLO-in-ag approach 
it is the small x - large y region where (^^jQ^) is beginning to systematically undershoot data. 
Hence, correcting the values of F2^{x, Q^) in the fit above for the correct values of Fl{x, Q^) leads 
to a worsening of the fit by ~ 25 (even after the fit is redone) and a consequent reduction of the 
confidence level to 0.5%. This effect, which is always ignored in LO-in-ttg fits to structure function 
data, clearly weakens the support for this type of approach considerably. 

It is when I perform the consistency check between the values of F|"(x,Q^) and F£(x,(5^), 
i.e. examine those points which are most sensitive to Fl{x,Q'^), that the most dramatic effects 
are seen. The values of -Fg (x, Q^) are fairly similar to those in the NLO-in-ag case, and hence so 
are the values of F^(x,Q^). Also, any significant change in the extrapolated Ff(x,Q^) reduces 
the confidence level of the fit to unacceptable levels, so once again the error on F^{x, Q"^) is taken 
to be that on a{x,Q'^) alone, flowever, the values of F^{x^Q'^) are very much larger than in the 
previous approaches. The set of values of F^{x,Q'^) and F^{x,Q'^) are shown in table 5 and 
displayed in fig. 4 . It is immediately clear that the measured values undershoot the predicted 
values significantly. A calculation of the is not necessary to prove that an extrapolation into 
the region where both F2{x, Q^) and Fl{x, Q"^) make important contributions to the cross-section 
rules out a LO-in-ag calculation of structure functions completely. 

An alternative way to demonstrate this is to compare the values of F|'(x,(5^) obtained from 
the fit to the low y data to the "measured" values F^(x,Q^) obtained from the measurement of 
the cross-section and the predicted values of F^(x, Q^). This is shown for two x bins in fig. 5, and 
it is very easy to see that the large predicted F^{x, Q^) combined with the measured a(x, Q"^) leads 
to a large upturn in the data, above the theoretical curve, as y becomes larger. A very similar 
effect is seen in all the different x bins, and is clear evidence that a LO-in-ag fit is completely ruled 
out by high y data. I note that in this LO-in-CKg fit a large component of the small x F2{x,Q^) 
comes from the steepness of the input quark. If I were to follow the reasoning behind double 
asymptotic scaling, i.e. that the small-x inputs should be flat at low Q^, and the rise generated by 
evolution, then not only would the quality of the global fit fall, but the value of F£(a;, Q^) would 
be even larger, and this clear discrepancy between theory and data at high y would become even 
more pronounced. Therefore, this approach is very strongly ruled out in the region of small x and 
large y. The fact that the LO-in-ag approach fails badly, but acquires very important small-x 
corrections at NLO-in-ag also implies that the large small-x corrections at higher orders in 
should be important, although, as already mentioned, there is no overwhelming evidence for this 
yet. 

The three different theoretical approaches used above are far from being the full set proposed 
to describe small x physics, and this technique can be applied to others. However, perhaps rather 
alarmingly they are the only approaches which have been used in fits to a global range of data and 
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which are fully constrained (except for approximate fits using factorization-scheme-dependent, 
and hence incorrect methods of including leading ln(l/x) terms within the collinear factorization 
framework), and considerable work would be required to use any other approach. 

I also note that at rather higher x than considered in this paper there have been a number 
of direct measurements of Fl{x,Q^): the SLAC hydrogen and deuterium scattering experiments 
[16], the SLAC E140 experiments [17], CDHSW [18], BCDMS [19], CCFR [20] and NMC [11]. This 
data has fairly large error bars, and much is in the region of the charm threshold. However, it 
seems that with a better treatment of this threshold than used in this paper then either the NLO- 
in-a^ approach or the LORSC approach would fit the data fairly well, with the former tending 
to be perhaps a little high, and the latter perhaps a little low. (The RqCD curve in [11] should 
be treated with caution since it uses the gluon from [1], which has been criticised above, and a 
LO-in-a^ formula using four massless quarks rather than the appropriate NLO-in-CKg formula with 
massive charm quarks.) The very low data from these measurements gives some evidence for 
higher twist effects [21]. 

I finally note that during the preparation of this article some new preliminary data at high 
y has been released by HI [22]. There is new data at y = 0.7 for the same values of Q^, except 
for the lowest bin, which seems to improve the comparison for both the NLO-in-ag fit and the 
LORSC fit, i.e. the measured cross-section for = 12GeV^ and 20GeV^ increase, while that for 
25GeV^ goes down a little. The improvement for the LORSC fit appears to be better than that 
for the NLO-in-Os fit. The comparison for the LO-in-CKg will get even worse. There is also some 
data at y = 0.82 with large error bars. Here the cross-section measurements look rather low. The 
QCD study is presumably performed in the same manner as in [4]. A detailed study in the manner 
of this paper will await final data. 

3. Conclusions. 

In this paper I have examined the implications of the idea proposed in [4] of extrapolating a 
fit performed to structure function data at relatively low y into the region of high y, where the 
longitudinal structure function starts to make an impact on the cross-section. In [4] it was assumed 
that believing the extrapolation of F2(x,Q^) into this high y region using a given theoretical 
framework, the difference between the extrapolation and the measured cross-section would give a 
measurement of F2{x,Q^). In this paper I have worked on a different principle, pointing out that 
different theoretical models will lead to different extrapolations, and also different predictions for 
Fl{x,Q^). All one is really examining is whether, in the region where both structure functions 
contribute to the cross-section, the data on the cross-section are consistent with the theory. One 
way to present this is to compare the "measured" values of Fl{x,Q'^) with the predicted values. 
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I have first done this using the same theoretical framework as in [4] (up to a different treatment 
of charm mass effects) and fitting to the same "low y" data. This provides results almost identical to 
[4], which may be naively interpreted as consistency of the NLO-in-a^ structure functions at high 
y. However, I have also argued very strongly that this type of fit is very badly under constrained, 
particularly concerning the gluon for a; > 0.1, and is not useful. I have then repeated the procedure 
using a NLO-in-a^ global fit at low y which is a great deal more constraining. I have also used 
points in the latter of [2], and also in [1], which give similar uncertainty on an extracted Fl{x, Q^) 
to those points in [4]. Even using this same theoretical framework, the conclusion concerning the 
measured values of Fl{x,Q'^) changes somewhat, with the average value being ^ 60% those in [4]. 
The agreement between the "measured" and predicted values is good. 

Using two further theoretical approaches, the LORSC calculation and the LO-in-a^ calculation 
fits to the same low y data the "measured" values of Fl{x, Q^) do not change by very large amounts, 
although the changes would quickly increase with increasing y. The predictions for Fi,{x,Q^) are 
however different to the previous case: for the LORSC approach consistency is in fact very slightly 
better than the NLO-in-a^ approach, but in the LO-in-ag approach the predicted Fl{x, Q^) is far 
too large. Thus, the LO— in-ctg approach, and consequently (and particularly) double asymptotic 
scaling, is ruled out precisely in the region where it is claimed most strongly to hold, i.e. high 
and low x. 

Hence, within this limited study I conclude that both the NLO-in-Os and the LORSC calcula- 
tions are consistent with high y cross-section data. Therefore, forgetting other possible theoretical 
models, the values of F]^[x, Q"^) are very likely to be similar to the predictions of these approaches 
when they arc constrained by a global fit to structure function data. By definition this means that 
the values arc similar to those "measured" using the two approaches, but somewhat lower than 
those presented in [4]. An increase in the precision of measurement of the high y cross-section, or 
an extension to slightly higher values of y would be very important in differentiating between the 
two theoretical approaches, and potentially any others. 

Alternatively, a direct measurement of Fl(x, Q^) would be even more useful. In this case there 
would not be any inbuilt uncertainty due to a particular (and always to some degree approximate) 
theoretical model, or to the measurement depending on how accurately, or correctly, unknown 
parameters are determined by the fit to low y data (even though this is probably small) . It would 
also eliminate the significant possibility that incorrectness in both F2(x,Q^) and Fl{x,Q'^) in a 
given theoretical approach could act to partially, or even largely cancel out and still lead to apparent 
consistency in the above approach. Hence I encourage strongly any attempt to measure Fl{x,Q'^) 
directly using any method as a way to help discriminate strongly between different theoretical 
approaches to calculating small x structure functions, and to obtain real data on a real physical 
quantity. 
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Table 1 

Data points with sufficiently high y and low enough error on a for the error on Fp{x,Q'^) to be 
less than 0.3. The values are denoted in GeV^. 



HI data 


ZEUS data 


Q' 




y 


a 


Aa 


Q' 


X- 10^ 


y 


a 


Act 


8.5 


1.35 


0.70 


1.165 


0.099 


6.5 


1.00 


0.72 


1.092 


0.084 


12.0 


1.90 


0.70 


1.198 


0.079 


8.5 


1.60 


0.59 


1.248 


0.074 


15.0 


2.38 


0.70 


1.368 


0.085 


10.0 


1.60 


0.69 


1.225 


0.114 


20.0 


3.17 


0.70 


1.276 


0.079 


12.0 


2.50 


0.53 


1.233 


0.068 


25.0 


3.96 


0.70 


1.439 


0.089 


15.0 


2.50 


0.67 


1.424 


0.108 


25.0 


5.00 


0.56 


1.430 


0.069 


27.0 


6.30 


0.48 


1.407 


0.053 


35.0 


5.54 


0.70 


1.435 


0.099 


35.0 


6.30 


0.62 


1.459 


0.086 












60.0 


10.0 


0.67 


1.438 


0.088 



Table 2 

The "measurement" Fj^{x,Q^) and prediction F^{x,Q^) of Fl{x,Q^) for the NLO-in-a^ fit per- 
formed in a manner analogous to that in [4], and for the data in [4] only. Also shown are the 
"measured" values F£'{x,Q'^)[4] in [4]. 



HI data 




X ■ 10^ 






AF}^ 


F-[4] 


8.5 


1.35 


0.31 


0.51 


0.22 


0.51 


12.0 


1.90 


0.34 


0.62 


0.18 


0.63 


15.0 


2.38 


0.35 


0.33 


0.19 


0.35 


20.0 


3.17 


0.36 


0.61 


0.17 


0.67 


25.0 


3.96 


0.39 


0.30 


0.20 


0.33 


35.0 


5.54 


0.38 


0.33 


0.22 


0.39 
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Table 3 

and Fl{x,Q'^) for the NLO-in-a« global fit. 



HI data 


ZEUS data 




a; -10^ 




Fr 




Q' 


X- 10^ 


Fl 


Fr 




8.5 


1.35 


0.24 


0.41 


0.22 


6.5 


1.00 


0.19 


0.39 


0.17 


12.0 


1.90 


0.28 


0.51 


0.18 


8.5 


1.60 


0.23 


0.16 


0.25 


15.0 


2.38 


0.30 


0.21 


0.19 


10.0 


1.60 


0.26 


0.34 


0.26 


20.0 


3.17 


0.31 


0.50 


0.17 


12.0 


2.50 


0.270 


0.42 


0.29 


25.0 


3.96 


0.35 


0.19 


0.20 


15.0 


2.50 


0.29 


0.03 


0.27 


25.0 


5.00 


0.33 


0.02 


0.27 


27.0 


6.30 


0.32 


-0.18 


0.30 


35.0 


5.54 


0.38 


0.22 


0.22 


35.0 


6.30 


0.33 


0.05 


0.26 












60.0 


10.0 


0.32 


0.07 


0.22 



Table 4 

F|"(x,Q2) and Fl{x,Q^) for the LORSC global fit. 



HI data 


ZEUS data 




X- 10'^ 


Fl 


Fr 


AF™ 




X ■ 10'^ 


Fl 


Fr 




8.5 


1.35 


0.20 


0.43 


0.22 


6.5 


1.00 


0.19 


0.37 


0.17 


12.0 


1.90 


0.22 


0.53 


0.18 


8.5 


1.60 


0.19 


0.12 


0.25 


15.0 


2.38 


0.22 


0.25 


0.19 


10.0 


1.60 


0.21 


0.33 


0.26 


20.0 


3.17 


0.22 


0.55 


0.17 


12.0 


2.50 


0.20 


0.37 


0.29 


25.0 


3.96 


0.24 


0.24 


0.20 


15.0 


2.50 


0.22 


0.03 


0.27 


25.0 


5.00 


0.23 


0.09 


0.27 


27.0 


6.30 


0.21 


-0.18 


0.30 


35.0 


5.54 


0.23 


0.30 


0.22 


35.0 


6.30 


0.22 


0.09 


0.26 












60.0 


10.0 


0.21 


0.13 


0.22 
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Table 5 

and Fl{x,Q'^) for the LO-in-a« global fit. 



HI data 


ZEUS data 




ar- 10^ 




Fr 




Q' 


10^ 


Fl 


Fr 




8.5 


1.35 


0.61 


0.35 


0.22 


6.5 


1.00 


0.60 


0.33 


0.17 


12.0 


1.90 


0.60 


0.45 


0.18 


8.5 


1.60 


0.58 


0.07 


0.25 


15.0 


2.38 


0.59 


0.17 


0.19 


10.0 


1.60 


0.61 


0.28 


0.26 


20.0 


3.17 


0.57 


0.47 


0.17 


12.0 


2.50 


0.56 


0.33 


0.29 


25.0 


3.96 


0.60 


0.17 


0.20 


15.0 


2.50 


0.58 


-0.01 


0.27 


25.0 


5.00 


0.56 


0.01 


0.27 


27.0 


6.30 


0.51 


-0.18 


0.30 


35.0 


5.54 


0.55 


0.23 


0.22 


35.0 


6.30 


0.53 


0.07 


0.26 












60.0 


10.0 


0.45 


0.13 


0.22 
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Figure Captions 



Fig. 1. Result of the measured values of F^{x,Q^) for the six points used in [4], along with the 
predicted value, F^{x, Q^), along a curve of constant y = 0.7 (the solid line), when using a 
NLO-in-Q!s fit analogous to that in [4]. Also shown are the measured values of Fl{x, Q^) 
in [4], F^{x, Q'^)[4], and the curve of predicted values in this paper (dashed curve). 

Fig. 2. Result of the measured values of -F™(x, Q^) for the 15 points satisfying the cut of 0.3 on 
the error, along with the predicted values, Fl{x,Q'^), for a NLO-in-ag global fit. In the 
top figure the points all have y ^ 0.7 and Fl{x, Q^) is along a curve of constant y = 0.7 
(the solid line). In the lower figure the points have y 0.55 and F£(x, Q^) is along a 
curve of constant y = 0.55. 

Fig. 3. Result of the measured values of Fj^{x,Q^) for the 15 points satisfying the cut of 0.3 on 
the error, along with the predicted values, Fl{x,Q'^), for a LORSC global fit. 

Fig. 4. Result of the measured values of Fj^{x,Q^) for the 15 points satisfying the cut of 0.3 on 
the error, along with the predicted values, Ff^{x,Q'^), for a LO-in-ag global fit. 

Fig. 5. The "measured" values of the structure function, F^{x,Q'^) (data points), compared to 
the theoretical values, F2{x,Q^), obtained from the best fit (solid line) in two different 
small X bins, when using the LO-in-ag calculation. 
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